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ABSTRACT 



We report new mid- infrared observations of the remarkable object IRAS 
15099—5856 using the space telescopes AKARI and Spitzer, which demonstrate 
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the presence of prominent crystalline silicate emission in this bright source. IRAS 
15099—5856 has a complex morphology with a bright central compact source 
(IRS1) surrounded by knots, spurs, and several extended (~ 4') arc-like fila- 
ments. The source is seen only at > 10 /im. The Spitzer mid-infrared (MIR) 
spectrum of IRS1 shows prominent emission features from Mg-rich crystalline 
silicates, strong [Ne II] 12.81 /im and several other faint ionic lines. We model 
the MIR spectrum as thermal emission from dust and compare with the Herbig 
Be star HD 100546 and the luminous blue variable R71, which show very similar 
MIR spectra. Molecular line observations reveal two molecular clouds around 
the source, but no associated dense molecular cores. We suggest that IRS1 is 
heated by UV radiation from the adjacent O star Muzzio 10 and that its crys- 
talline silicates most likely originated in a mass outflow from the progenitor of 
the supernova remnant (SNR) MSH 15—52. IRS1, which is embedded in the 
SNR, could have been shielded from the SN blast wave if the progenitor was in 
a close binary system with Muzzio 10. If MSH 15—5,2 is a remnant of Type Ib/c 
supernova (SN Ib/c), as has been previously proposed, this would confirm the 
binary model for SN Ib/c. IRS1 and the associated structures may be the relics 
of massive star death, as shaped by the supernova explosion, the pulsar wind and 
the intense ionizing radiation of the embedded O star. 

Subject headings: infrared: stars — stars: individual (IRAS 15099—5856) - 
ISM: individual (MSH 15— 5£) — supernova remnants — circumstellar matter 



Introduction 



Space infrared observations have revealed distinct spectral features due to crystalline 



silicates in diverse objects such as Asymptotic Giant Branch (AG B) stars, p 



Herbig Ae/Be stars, comets and ultraluminous infrared galaxies (IHenning 



anetary nebulae, 



2010 ). Since sil i- 



cate dust grains in the interstellar medium are essentially amorphous (jKemper et al.l 12004 ). 
crystalline silicates must form in circumstellar disks and/or outflows of evolved stars or young 
stellar objects (YSOs). However, neither the formation process of the crystalline dust nor 
its relation to the central stellar source is understood. 

In this paper, we report the discovery of prominent crystalline silicate emission in IRAS 
15099—5856, a bright, mid-infrared compact source previously detected by Infrared Astro- 
nomical Satellite (IRAS) and Midcourse Space Experiment (MSX). The source is located close 
to the pulsar B1509-58 in the supernova remnant (SNR) MSH 15-53 (G320.4-1.2), a young 



SNR of complex morphology at a distance of 5.2 ±1.4 kpc (lArendtlll99ll ; lGaensler et al 



1999 
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see § 5). lArendtl (Il99ll ) concluded that IRAS 15099—5856 is heated either by hot plasma 
associated with th e SNR or by the n earby O star Muzzio 10 (jMuzziolll979l ). which is at a 
distance ~ 4 kpc (IBessell et al.ll201ll ). We show that the latter explanation is most likely 
correct and that the source is probably associated with the progenitor of the SNR. 



2. AKARI Images and Spectral Energy Distribution 

We observed IRAS 15099—5856 using the Infrared Camera (IRC) aboard AKARI on 
February 21-22, 2007 (ID: 1400761.1 and 1400762.1). The IRC is equipped with three wave- 
band channels cove ring the 2.6-26 /im wavelength range with a 10' x 10' field-of-view (FOV) 



(jOnaka et al.l 120071 ). We used all three channels in IRC02 mode, which gave six band im- 
ages, N3, N4, S7, Sll, L15, and L24, with central wavelengths of 3.2, 4.1, 7, 11, 15 and 24 
fim, respectively, and angular resolutions ranging from 4. "2 (N3) to 6. "8 (L24). The data 
were processed using the standard IRC Imaging Data Reduction Pipeline version 070104. 
The positional uncertainty (la) is < 0."15 for N3- and N4-band images and < 0."32 for the 
longer-wavelength images except the L24-band image, which has a relatively large uncer- 
tainty (~ 1") because the astrometric solution could not be obtained by matching with the 
coordinates of 2MASS stars in the Pipeline. 

The AKARI observations reveal that IRAS 15099—5856 is a large, MIR source with 
spectacular morphology (Fig. 1). It is composed of a bright central compact source (hereafter 
IRS1), a surrounding halo of ~ V radius with knots and spurs embedded in it, and several 
extended (~ 4'), knotty arc-like filaments. The extent of the whole structure is about 10' 
and its total flux is ~ 12 Jy at 15 /im, 40% of which is from IRS1. The 15 /im image in 
the right frame shows that IRS1 has a brightest central part extended along the east-west 
direction and some diffuse emission around it. There is no 2MASS counterpart to IRS1. The 
source is not seen in the IRC images at wavelengths A < 10 /im. 

Figure 2 shows the one- dimensional 24, 15, and 11 fim intensity profiles of IRS1 along the 
east- west (P. A. =110°) and north-south directions, which are the two orthogonal directions 
of the AKARI IRC. We have subtracted a constant background brightness from each profile 
and normalized it by its maximum brightness: background brightness of about 150, 60, and 
30 MJy sr _1 and maximum brightness of 3700, 1110, and 30 MJy sr _1 for L24, L15, and 
Sll profiles, respectively. For comparison, we also show the point spread function (PSF) 
of the AKARI IRC. The 15 and 24 /im profiles clearly show that IRS1 is composed of a 
compact bright central part (hereafter, the core) and an extended envelope of ~ 30" size. 
The envelope appears as a plateau in the north-south profile, which suggests that it is a 
separate component from the core. The core is almost symmetrically elongated along the 
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east-west direction, while it is barely resolved in the north-south direction. The average 
full- width at half- maximum (FWHM) of the core at 15/24 /im after deconvolving the PSF 
is 9. "6 (0".l) x 5."1 (l."l) or 0.19^4 pc x0.10d 4 pc, where d 4 is the distance in units of 4 kpc. 
(We normalize the distance by 4 kpc, which is the distance to the Cir OBI association in 
this area. See § 5.) At 11 /im, the core has a comparable size (9. "7) along the east-west 
direction, but its center is slightly (l."5) shifted to the west. Along the north-south direction, 
it appears somewhat more extended (6. "4), but this could be due to the contribution from 
nearby faint stellar sources (see Fig. 1). The local peak at —13" in the 11 /im north-south 
profile is due to Muzzio 10. 

Table 1 summarizes the observed fluxes of IRAS 15099—5856, where we also list the 
corresponding "color-corrected" fluxes, i.e., the fluxes corrected for the spectral slope over the 
passbands. The color correction has been made by assuming modified black-body emission 
at 75 K that fits the overall shape of the spectral energy distribution (SED). But as we show 
in § 3, the MIR emission at < 35 /im is dominated by spectral bumps from crystalline and 
metallic dusts, so that the color-corrected fluxes of those bands should be used with caution. 
The AKARI IRC fluxes are derived by summing the fluxes within a circle of 25" radius 
centered on the brightness peak; the contribution from the O star Muzzio 10 is excluded. 
The IRAS and MSX fluxes are from the point source catalogs and they roughly agree with 
the AKARI IRC fluxes. At wavelengths shorter than the AKARI IRC bands, we identified 
several point sources of 17 th -20 th mag around IRS1 in the K s image of the source obtained 



with the Magellan 6.5 m telescope (IKaplan &: Moonll2006l ). They are separated from the 15 



/xm peak position by > I'M, which is greater than the 3a positional uncertainty of the 15 /im 
image. We use the 3a upper limit (20 th mag or 0.007 mJy) of the Magellan observations as 
an upper limit at 2.15 /im. At longer wavelengths, we identified the source in the AKARI 
FIS (Far-Infrared Surveyor) catalog at 65 and 90 /im but not at 140 or 160 /im. With a 
spatial resolution of 37"-61", the FIS does not resolve IRS1. We also carried out 1.1 mm 
continuum observations using the 144-element bolometer array AzTEC ( Wilson et al. 2008h 



on th e 10-m telescope Atacama Submillimeter Telescope Experiment (ASTE ; 



Ezawa et al. 



20041 ) in October 2008. The beam size and field of view were 28" and 7/5, respectively. 
We obtained 11 maps of a 25' x 25' area in a raster scan at 90" per sec, and coadded the 
individual maps to obtain the final map with a pixel size of 6". The total on-source observing 
time was 3 hrs and the per-pixel integration time was 10 s. The absolute flux calibration 
was done by observing Uranus. IRAS 15099—5856 was not detected with an upper limit of 
30 mJy (2a). 

Figure 3 shows the SED of IRS1. The flux of IRS1 increases steeply from short wave- 
lengths to 25 /im, and remains flat at ~ 30 Jy to 90 /im. The MIR (11 — 25 /im) spectral 
index is a = —7.7 (F v oc v a ). The SED can be approximately fitted by emission from 
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2.4 x 10~ 3 d1 M & of dust at 79 K using Dra ine's carbonaceous-silicate opacity model for 
interstellar dust with Ry = 5.5 (jDraine 2003 ). The fit, however, is not good, with a large 
reduced x 2 {= 9.4). This is because the flux at £j 35 /im is mainly from dust of circumstellar 
origin with several distinct opacities (§ 3). In § 3, we discuss the model based on the Spitzer 
IRS spectra, where we show that the dust mass is about four times larger because the dust 
temperature is lower. The total IR luminosity of IRS1 is Lirsi ~ 2.0 x lO 3 ^ L Q . 



3. Spitzer IRS Spectrum and Crystalline Silicate Emission 



We carried out Spitzer IRS staring observations of IRS1 of IRAS 15099—5856 with the 
two low-resolution modules; the Short-Low (SL) module (5.2-14.5 /im) and the Long-Low 
(LL) module (14.0-38.0 /im). The observations were done on October 3, 2008 (ID 26318080 
and 26318336). We followed the standard IRS staring mode, where the spectra were obtained 
by placing the source at two locations for each slit, located 1/3 and 2/3 of the way along 
the slits' length (Fig. 1). At each location, 2-3 cycles of different ramp durations (6 s to 
60 s) were carried out. We use the coadded 2-D spectra provided by the Spitzer Science 
Center (SSC), which are averaged basic calibration data (BCD) FITS files processed by 
the pipeline version 17. For the data analysis, we use the Spitzer IRS Custom Extraction 
(SPICE) software v2.2 provided by the SSC. The IRS1 is extended, so that the flux that 
misses the slit, i.e. "slit-loss", needs to be properly accounted for the flux calibration (see 
the SSC IRS calibration web page on extended sources^). We derive the slit-loss correction 
factor (SLCF) at 11, 15, and 24 /im by calculating the fraction of t he flux with i n the S L and 
LL slits on the AKARI images using the Point Spread Function in lSloan et al.l (j2003a[ ). The 
SLCFs over the entire SL and LL bands are then obtained by interpolating/extrapolating 
these data points. 

The Spitzer IRS spectrum of IRS1 (Fig. 4) increases steeply from £ 0.1 Jy at £ 13 /im 
to 20-30 Jy at > 20 /im, which is consistent with the broadband SED. The spectrum shows 
two broad, strong peaks at 23 and 34 /im and also weak features at 20 and 27 /im. Such 
spectral features are characteristic of Mg-rich crystalline silica tes, e.g., forster ite (Mg 2 Si04), 
and have been observed in AGB stars, YSOs, and comets (iHenningl 120101 ). Indeed, the 
> 15 /im-part of the spectrum is very similar to those of YSOs such as HD 100546 or 
RECX 5 (Fig. 4). HD 100546 is the Herbig Be star (B9.5) that has an exceptionally 
l arge crystallinity in the Infrared Spa ce Observatory (ISO) sample of Herbig Ae/Be stars 
f lMalfait et aTlll998l : iMeeus et al.lboOlh . RECX 5 is an M4-type pre-main-sequence star in the 



http://ssc.spitzer.caltech.edu/irs/calib/extended_sources/ 
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?7 Ch amaeleontis cluster and shows a very similar spectrum to HD 100546 (jBouwman et al. 
20051 ). The spectrum is also quite similar to that of luminous blue variables (LBVs) such 
as R71 (HD 269009) in the Large Magellanic Clou d (LMC) (Fig. 4) . R71 is a well-studied 
LBV with T* = 15,000 K and L = 5 x 10 5 L faover et al.lboioh . Its ISO and Spitzer 
IRS spectra show a broad bump ste eply rising from 15 /im and flattenin g; around 18 jum, 
which is very similar to that of IRS1 ( Voors et al. 1999[ Morris et al. 2008 ). This bump had 
been attributed to the metal oxide FeO (IVoors et al.lll999l ). It also shows the weak bumps 
due to Mg-rich crystalline silicates. The closely-matching spectral features suggest that the 
grain composition should be similar in these sources. A noticeable difference, however, is the 
absence of crystalline/amorphous silicate emission at 10-11 jum in IRAS 15099—5856, which 
is distinct in the other sources. This indicates that IRAS 15099—5856 is missing the warm 
(several hundred K) silicate dust that the other sources have. Before closing the comparison 
to other MIR spectra with similar features, we note that dust emission has been recently 
detecte d toward several Oxygen-rich SNRs, e.g Cas A, G292.0+ 1.8, and the LMC SNR 
N132D faho et alil2"o~0~8l iGhavamian et al.lEooalTappe et al.ll2006h . The 17-22 /im features 
in the spectra of these objects have been attributed to either silicate dust newly formed in 
supernova ejecta or swept-up polycyclic aromatic hydrocarbons. They, however, do not show 
crystalline silicate features. 

Several ionic lines are present in IRAS 15099—5856: strong [Ne II] 12.81 jum and weak 
[Ne III] 15.56 /im, [S III] 18.71 /im, and [O IV] 25.89/[Fe II] 25.99 /im. The central velocity 
of the [Ne II] line is -160 ± 560 km s _1 , which indicates that IRAS 15099-5856 is not at 
a cosmological distance. The flux of the [Ne II] 12.81 /im line is ~ 1.2 x 10~ n erg cm -2 
s _1 . The ratio [Ne III]/[Ne II] ~ 0.14, although this is uncertain since the lines appear in 
different, almost orthogonal, slits. The spatial extent (FWHM) of the [Ne II] line-emitting 
region along the north-south-oriented SL slit is 5", while that of the dust-continuum-emitting 
region along the east-west oriented LL slit is 12"-13". Considering the uncertainties in the 
size measurement with the slit spectrograph, e.g., slight misalignment of the slit from the 
peak position, we consider that these results are consistent with the size of IRS1 determined 
from the AKARI IRC image in § 2. 

We model the IRS spectrum as thermal emission from several independent dust compo- 
nents. We do not carry out a detailed analysis considering all possible dust species, which is 
beyond the scope of this paper. We instead determine the minimum number of dust species 
required to explain the observed IRS spectra and derive their characteristic parameters. 
The IRS spectrum in Figure 4 requires at least three dust components. First, the narrow 
intensity peaks at 22.6, 27.2, and 33.5 /im clearly indicate a considerable amount of Mg-rich 
crystalline olivine, i.e., Mg2 X Fe2-2xSi04 with x ~ 1. The peak frequency and the relative 
strength of the mass absorption coefficient (hereafter MAC) depend on the Fe/Mg ratio, 
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Fabian et al. 


2001; 


Koike et al. 


2003, 


2006) 



The strong p eaks at the above wavelengths are characteristics of Mg-rich, Fe-poor olivines 
(see Fig. 1 o f Koike et al.ll2003l ). We adopt Mgx 9 Feo.iSi0 4 and use the optical constants of 



Fabian et al.l (j200lf ). We assume spherical dust grain s with uniform radii of 0.1 /j,m and calcu- 



late their MAC following iBohren fc Huffman! (119831 ). Second, the broad bump steeply rising 
from 15 fim and flattening around 18 /im could be attributed to metal oxides (Fe x Mgi_ x .O, 
< x < 1). Metal oxides have single-peaked, broad MACs with p eak wavelengths between 



16 and 22 /im depending on the Fe content (IHenning et al. 



1995). We use a combination 



of FeO and Mgo.6Feo.4O with optical constants of IHenning et al.l (119951 ). The MACs of 0.1 
/im spherical particles, however, are rather narrow to fit the observed broad bump, and we 
use the MAC for continuous distribution of ellipsoids for FeO using a quadratic weighting 
where near-spherical particle shapes are most probable. Different combinations of metal 
oxides might fit the spectrum equally well. For example, MgO has an absorption peak at 
wavelengths shorter than that of FeO and can be used instead of Mgo.6Feo.4O (see below). 
Finally, we need an additional, cool dust component for the far-IR emission that peaks at 
about 60 /im. Since we have only broad-band fluxes at > 35 /im, we cannot constrain the 
composition of this component. We adopt a morphous silicate, and use the optical constants 
of cool silicate grains derived by ISuhl (1l999h . 



With the above four dust components, i.e., crystalline olivine, FeO, Mgo.6Feo.4O, and 
amorphous silicate, we fit the Spitzer IRS spectrum using MPFIT, which is a least-squares 
fitting tool based on the Levenberg-Marquardt algorithm (jMarkwardt 2009 ). In the fit, 
we leave the temperatures and the masses of individual components free with reasonable 
lower and upper boundaries. In order to constrain the far-IR part of the spectrum, we 
use the AKARI 65 and 90 /im fluxes. The best-fit masses and temperatures obtained 
by least-squares fitting are given in Table 2, where we also show the results when MgO is 
used instead of Mgo.6Feo.4O (Model 2). The resulting model spectrum fits the observed IRS 
spectrum and also the far-IR broad-band fluxes reasonably well (Figs. 2 and 5). The other 
model (Model 2) can fit the spectrum equally well (Fig. 5). We could not exactly reproduce 
the observed spectrum because of the simplicity of the adopted dust model and because the 
dust components are unlikely to be isothermal. The total dust mass is 9 x lO -3 ^ M Q , which 
is considerably greater than the estimate in § 2 because of the lower temperature of the 
silicate dust components. Most of the mass comes from amorphous silicate, the cool dust 
component for the far-IR emission, while crystalline silicates contribute about 25%. The 
temperatures and therefore the masses of these two dust components vary little between the 
two models. However, the parameters of the metal oxide components vary significantly, as 
can be seen from Table 2, so they are not reliable. T he temperature of c rystalline silicate (55 
K) is much lower than that (210 K) in HD 100546 iMalfait et alill998h . The corresponding 
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temperature in R71 is not available, but the temperature of the dust r esponsible for mos t 
of the emission at 10-100 /im in R71 has been estimated to be 120 K (IBoyer et al.ll2010l ). 
In R71, it w as noted that an a dditional warm dust at 300 K is needed to explain the bump 
near 10 /im (IVoors et al.lll999l ). 



4. Molecular Emission Lines and Other Wavebands 

We carried out molecular line observations of IRAS 15099—5856 using the 22-m telescope 
at the Mopra Observatory and ASTE . With Mopra, we mapped a 6' x 6' area surrounding 
IRAS 15099-5856 in the CO J=l-0 line at 115.27 GHz (FWHM=33") on June 28 and July 
13, 2008. With ASTE, we made 5-point cross observations at 20" spacing in the CO J=3-2 
line at 345.79 (FWHM=22") from 31 May to 6 June 2008. We detected molecular gas at two 
LSR velocities, —62 km s -1 and +22 km s -1 surrounding IRAS 15099—5856 (Figs. 6 and 
7). Neither component peaks at the position of IRS1. Instead, the —62 km s _1 -component 
shows a cavity at the position of IRS1 while the +22 km s _1 -component shows no obvious 
spatial correlation. The two clouds have a CO (J = 3 — 2)/(J = 1—0) ratio ~ 0.5, which 
is typical for giant molecular clouds. 

We have also searched the 1.4 GHz continuum data of this area obtained from the 
Australian Compact Telescope Array (Gaensler et al. 1999) and a deep (250 ks) Chandra 
X-ray archival observation, but could not identify a counterpart. The upper limit (5<r) on 
the 1.4 GHz continuum flux is ~ 10 mJy. 



5. Nature of IRAS 15099-5856 

The strong crystalline silicate emission indicates that the dust in IRS1 is circumstellar 
or perhaps of supernova origin. The steeply-rising MIR spectrum, the low-velocity [Ne II] 
12.8 /im line, and the large extinction through the Galactic plane in this direction almost 
rule out an extragalactic origin for IRAS 15099—5856. The derived dust mass (Table 2) and 
the spatial extent (0.19 x 0.10g? 4 pc) imply a small self-extinction ( ;S 0.1 mag at K band); 
together with the absence of an embedded star in images of IRAS 15099—5856, this rules out 
the possibility that it is an extended dust clump powered by an embedded star. It is unlikely 
that several highly-obscured stellar sources form IRS1, since the temperature (55-57 K) of 
the major dust components is less than expected for such circumstellar dust. Instead, the 
[Ne II] and [Ne III] emission lines and the IR luminosity suggest that it could be powered by a 
nearby hot, luminous star. According to stellar atmosphere models, the ratio [Ne III]/[Ne II] 
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~ 0.14 implies a stellar temperature of 30,000-40,000 K (IMorisset et al.ll2004j ). The observed 
[Ne II] luminosity i mplies a Lyman con tinuum photon luminosity Sirsi ~ 1 v 1 "~ J 
incident upon IRS1 f lHo and Ketol 120071 ). 



1.2 x 10 47 ^ s 



The nearby O star Muzzio 10 can provide both the req uired luminosity and the ionizing 
photons. Based on new optical spectra, iBessell et al.l ( 120111 ) concluded that its spectral type 
is 04.5III(fp) and it is probably a member of the Cir OBI association at ~ 4 kpc. We adopt 
T* = 40, 500 K, L, = 5.8 x 10 5 a nd = 3.2 x 10 49 s" 1 for Muzzio 10, the parameters of 
an 04.5111 star (jMartins et al.l 120051 ) . Note that the stellar temperature is consistent with 
that needed to provide the observed [Ne III]/[Ne II] line ratio. The projected distance of 
IRS1 from Muzzio 10 is 13". 7, or 0.27^4 pc, comparable to the size of IRS1. At 4 kpc, 
LrRsi/£* = 3.4 x 10~ 3 <i| and Sirsi/S 1 * = 3.8 x lO -3 ^, which are in excellent agreement. If 
the projected surface area of IRS1 as seen from Muzzio 10 is the same as that seen from 
Earth, then they are separated by 0.6 pc. The e xpected temperature of amorphous silicate 
grains of radius 0.1 /im is 70 K ( iTielensI 120051 ). but the mean temperature of dust in a 
cloud with a large UV optical depth will be lower. The temperature of 57 K inferred in 
§ 3 (Table 2) therefore is consistent with our expectation. For metal oxide dust, however, 
the expected temperature is considerably lower than the derived temperature (150 K) in 
§ 3. This, however, could be due to our assumed dust composition, e.g., MgO instead of 
Mgo.6Feo.4O yields 90-120 K instead of 150 K (Table 2). We conclude that the agreement 
between the predicted and observed Ne line luminosities, IR luminosities and silicate dust 
temperatures supports the hypothesis that IRS1 is powered by Muzzio 10, although the 
high temperature of the small amount of the metal oxide dust remains unexplained. This 
association implies that IRS1 cannot b e too old, since the ioniz ing radiation from Muzzio 10 
would have destroyed it in ~ 10 5 yrs (IBertoldi fc McKedll990l ); it would also be subject to 
destruction by the stellar wind from Muzzio 10. 

What is the origin of the crystalline silicate? It can be produced in both Y SOs and 
evolv ed stars, since it requires temperatures ^ 1000 K to form by thermal annealing ( IHenning 
20101 ) . The dust in YSOs is confined to circumstellar disks and is thus very compact; hence, 
a YSO origin does not apply to the present source, which is extended. The absence of 
dense molecular material associated with IRAS 15099—5856 is consistent with this conclu- 
sion. There are no luminous evolved stars within the boundaries of IRS1: Muzzio 10 is 
very hot and must still be close to the main sequence; it is not in the blue-loop phase of 
late-stage evolution of massive stars, and its location in the H-R diag ram is far from the 
"S Dor instability strip" that defines the left boundary of the LBVs ( iVinkl 120091 ) . There 
is a Wolf Rayet star (WR65; Oskinova & Hamann 2008) and another possible mass- losing 
star in the field, but they are relatively distant (4/3 south of and 4'.0 north of IRS1 respec- 
tively) and show very different MIR colors (see Fig. 1). The only candidate for the source 
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is the SNR MSH 15—5,2, which is at a comparable distance (5.2 ± 1.4 kpc — iGaensler et al. 
19991 ). or its progenitor. If we fit the broadband SED of Muzzio 10 from optical to MIR 
assuming a spectral type of 04.5111, we obtain a distance of 4.5 kpc and an absorbing H- 
nucleus column density of 9.2 x 10 21 cm -2 . This column de nsity matches well w ith that 



to the SNR determined from X-rays, 9.5 ± 0.3 x 10 21 cm -2 (IGaensler et al.l 120021 ) . which 
supports the association. The crystalline silicates could have been produced either in the 
late stages of the progenitor's evolution or in the cooling SN ejecta. We favor a progenitor 
origin for IRS1 in view of its low velocity, —160 ± 560 km s _1 . An origin in the progenitor 
is consistent with the ubiquity of crystalline silicates in evolved stars. Crystalline silicate 
emission features have not been detected in the spectra of freshly- formed SN dust, although 
crystallin e silicates of core-collapse SN origin have been discovered in interplanetary dust 
particles ([Messenger. Keller, fc Lauretta! 120051 ). Crystallin e silicates are sub ject to destruc- 
tion and amorphization by SNR shocks and cosmic rays (IJager et al.ll2003l . and references 
therein), so that they will not survive long in the hostile environment of a young SNR even 
if they are formed in the SN ejecta. 

A progenitor origin for IRS1 has several interesting consequences. First, IRS1 should 
be embedded in the SNR and be relatively cl ose to the explosion center, because otherwise, 
considering the large radius (23 pc at 4 kpc; IGaensler et al.lll999l ) of the remnant, the dust 
mass lost by the progenitor would be unacceptably large. This raises the question of how 
IRS1 could have survived the S N blast wave, because most dust g rains are destroyed by 
shocks faster than ~ 200 km s _1 (I Jones. Tielens. fc Hollenbachlll996l ). One possibility would 
be that the explosion was very asymmetric so that IRS1 avoided the expanding dense ejecta. 
Another, more intriguing possibility is that the pro g enitor and Muzzio 10 formed a close 
binary system, as first proposed by IGaensler et al.l (Il999l ). and that Muzzio 10 shielded 
IRS1 from the SN blast wave as shown, for example, in numerical simulations of Type la 
SN explosions in binary systems ([Marietta et al.l 120001 ) . The close proximity of Muzzio 10 
and the pulsar B1509— 5 8 gives support to t his possibility. Since Muzzio 10 is currently 
an 04.5111 star of 45 M ? (jMartins et al.l 120051 ) . the progenitor should have had a very large 
(60-70 M & ) initial mass and a separation ~ 100i? Q from Muzzio 10 (S.-C. Yoon, private 
communication). The progenitor could have lost most of its mass during Wolf-Rayet phase 
prior to its explosion as an SN. The shielding scenario for IRS1 implies a distance from the 
SN: If we assume that the diam eter of Muzzio 10 in a binary was comparable to that of an 
isolated 04.5111 star, i.e., 30i? Q (IMartins et al.ll2005l ). then it would shield a region 0.2 pc in 
extent (the observed size of IRS1) from a supernova 100-R© away if IRS1 were about 0.6 pc 
from the SN. This is the same as the estimate of the current separation between Muzzio 10 
and IRS1 made above. Dus t at thi s dista nce from the supernova could readily survive the 
SN, as can be inferred from iDraind (120101 ). If a significant fraction of the mass loss that led 
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to the formation of IRS1 were in a disk of that radius, then the total dust mass of the disk 
would be about (27r x 0.6 pc)/(0.2 pc) times the mass of IRS1, or about 0.17M Q . For a gas 
to dust ratio of 100, this implies a total mass of 17M , which is large but not unreasonable. 

Measurement of the proper motions of the pulsar and Muzzio 10 would be the key 
test of this binary shielding model. The orbits of the neutron star and its companion after 
the supernova explosion depend on their separation, masses, and the magnitude of the kick 
velocity imparted to the neutron star during the explosion (ITauris fc Takensl 119981 ). Since 
Muzzio 10 is much more massive than the neutron star, its velocity would be almost un- 
changed by the SN, and as a result its motion should have been nearly perpendicular to the 
vector from the SN to IRS1. The motion of the pulsar, on the other hand, could be mainly 
determined by the magnitude of the kick velocity and so could be in any direction. The 
available radio data on the pulsar B 1509— 58, however, do not show an appreciable proper 
motion (M. Livingstone, personal communication). Determinati on that the progenito r was 
a member of a binary system, together with the suggestion by iGaensler et al.l ( 119991 ) that 
the supernova w as of Type Ib/c, would also support the binary model for SN Ib/c (e.g., 
Yoon et al.l 120101 ). although this would require revision of the canonical vie w that only fas t 
winds are expected in the final evolutionary stages of Type Ib/c progenitors ( )Chevalierll2005l ) . 
Above, we commented that the mid-IR spectrum of IRAS 15099—5856 was similar to that 
of an LBV, which e xperiences e xtreme ma ss loss before i ts exp losion and is associated with 
SN Iln supernovae ( Smith 2010 ). However, Bessell et al. ( 201ll ) has obtained optical spectra 
toward IRS1 and detected no Ha emission. The absence of emission from hydrogen in the 
mass lost from the progenitor supports a Type Ib/c origin rather than a Type Iln origin for 
the supernova. 

The nature of the diffuse, extended structures surrounding IRS1 is u ncertain. They 



are n ot correlated with the X-ray knots or arcs of the pulsar wind nebula (IGaensler et al. 



20021 ). Spitzer IRS spectra of these structures indicate that they are a mixture of dust- 
continuum dominated material like IRS1 and line-dominated, ionized gas. It is likely that 
the red/yellow knots, spurs, and filaments in Figure 1 are dominated by dust continuum, 
while the diffuse green structures are dominated by line emission, in which case the former 
might be the progenitor's circumstellar material while the latter are SN ejecta. In one faint, 
[Ne II]-emitting gaseous clump, we measured a central velocity ~ +1, 000 km s _1 , consistent 
with SN ejecta. The SL and LL slits, which are almost orthogonal to each other, provide only 
limited spectral and spatial information to reach further conclusions. The spatial relation 
of IRS1 and the surrounding structures to the SNR as well as to Muzzio 10 needs further 
study. 



IRAS 15099—5856 appears to be the first case in which crystalline silicates have been 
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observed to be associated with a SNR or its progenitor. IRS1 and the associated structures 
may be the relics of massive star death in a close binary system, as shaped by the supernova, 
the pulsar wind and the intense ionizing radiation of the embedded O star. If confirmed 
by observations of the proper motion of Muzzio 10 and pulsar B1509-58, we have a unique 
opportunity to study the interior of a young SNR containing crystalline dust and illuminated 
by a central ionizing source. 
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Fig. 1.— (left) Three-color image of IRAS 15099-5856 produced from AKARI Sll (blue), 
L15 (green), and L24 (red) images. North is up and east is to the left. The brightness is 
adjusted to show the faint emission features, which saturated the central part. The cross 
marks the peak position of the central compact source (IRS1) at 15 /xm, (a 2 oocb &2000) = 
(15 h 13 m 56 s . 32, — 59 o 07'40."9). Th e plus and diamond signs mark the positions of the pulsar 
B1509-58 flGaensler et aD l2002h and the O star Muzzio 10 (2MASS 15135520-5907516), 
respectively. The bright star embedded in the diffuse blue emission in the south is the 
Wolf Rayet star WR65 while the other bright star of similar morphology in the north is not 
previously known, (right) Magnified view of the central area in L15 band (upper panel) and 
the corresponding 2MASS image (lower panel). The Spitzer IRS slit positions are shown by 
rectangular boxes in the upper panel; large (148" x 10.7") solid one representing LL band 
slit position and two superposed small (50.4" x 3.7") solid ones representing the two nodded 
SL band slit positions. 
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Fig. 2. — One-dimensional, normalized intensity profiles of the central compact source IRS1 
at 24, 15, and 11 /im (solid lines). The left panels show the profiles along the east-west 
direction (P.A. = 110°) while the right panels show the profiles perpendicular to it. The (0, 0) 
position represents the peak position at 15 /im, and the filled circles represent the intensities 
at individual pixels. The 24 fj,m image has a relatively large uncertainty in astrometry, and 
therefore the profile is shifted to match the peak position at 15 /im. The pixel sizes are 
2. "34-2. "38. The dotted lines show the profiles of a point source in the same images. 
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Fig. 3. — Spectral energy distribution of IRS1 of IRAS 15099—5856. The open symbols 
represent upper limits. The fluxes are color corrected (see Table 1). The dotted line is a 
model SED of thermal dust emission at 79 K, whereas the solid line is a model SED resulting 
from a fit to the Spitzer IRS spectrum (Model 1 in Table 2; see § 3). 
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Fig. 4. — Spitzer IRS spectrum of IRAS 15099—5856. Wavelengths of the identified ionic 
lines are marked. The > 15 //m-part of the IRAS 15099 —5856 spectrum is very similar to 
the spectrum of the Herbig Be (B9.5e) star HP 100546 JMalfait et[all ll998[) and the LBV 
R71 (HD 269006) in the LMC JVoors et alill999l : iMorris et allboosh . which are shown for 
comparison. Th e HD 10054 6 spect rum is from the ISO Short Wavelength Spectrometer 
spectral atlas of ISloan et al.l (I2003bl ) while the R71 psectrum is from the Spitzer Heritage 
Archive. The bottom panel shows the relative responses of the AKARI S7, Sll, L15, and 
L24 imaging bands. 
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Fig. 5.— Model fits to the Spitzer IRS spectrum of IRAS 15099-5856 (thick solid lines); (a) 
Model 1 and (b) Model 2 in Table 2, respectively. The thin solid line in each plot represents 
a model of thermal emission from several independent dust components, which are shown. 
The filled-star symbols mark the 15 and 24 /zm AKARI fluxes. The bottom frames show 
the differences between the model fits and the observed spectrum. 
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Fig. 6. — CO J = 1 — and J = 3 — 2 line profiles toward IRS1. The former is obtained 
by the 22-m Mopra telescope (FWHM=33") while the latter is obtained by the 10-m ASTE 
telescope (FWHM=22"). 
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Fig. 7. — CO J = 1 — line intensity distribution of two molecular clouds detected to- 
ward IRAS 15099-5856; the -62 km s" 1 (-65 to -61 km s" 1 ) and +22 km s" 1 (+16 to 
+25 km s -1 ) clouds are in the left and right frames, respectively. Contour levels repre- 
senting integrated intensity are 20, 30, 40, 50, 70, 90, and 110 K km s _1 . The gray-scale 
plot distinguishes hills from valleys. The red contours show the 24 /zm brightness of IRAS 
15099-5856. 
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Table 1. Flux of IRAS 15099-5856 IRS1 



Flux 



Wavelength Observed Color-Corrected 



Telescope 


(/im) 


(Jy) 


(Jy) 


IVTaffellan 

J. V J. CAi CL Vj i A CAi X J. 


2.15 


< 7 x 10~ 6 




MSX 


8.28 


< 0.101° 






14.65 


1 43 (0 W) 


1 20 (0 08) 




21.3 


12.88 (0.77) 


12.50 (0.75) 


AKARI IRC 


7 


< 0.0013 






11 


0.116 (0.003) 


0.076 (0.002) 




15 


5.14 (0.15) 


1.75 (0.02) 




24 


21.4 (1.0) 


23.5 (1.1) 


AKARI FIS 


65 


27.9 (4.8) 


26.6 (4.6) 




90 


17.2 (2.5) 


11.0 (1.6) 




140 


< ll 6 






160 


< 25 6 




IRAS 


12 


< 1.89 






25 


19.9 (2.0) 


23.8 (2.4) 




60 


37.0 (4.1) 


31.5 (3.5) 




100 


< 138 c 




AzTEC 


1100 


< 0.03 





Note. — Magellan, AKARI IRC, and AzTEC fluxes are from 
this work. For the non-detections, the listed fluxes are 2a upper 
limits. The other fluxes are from the Point Source Catalogs. 
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a MSX flux limit with quality _flag=l flEgan et all 120031 ). 

b Detection limits according to the FIS Bright Source Catalog 
Version 1. 

c Upper limits from the IRAS Point Source Catalog Version 
2.0, which are nominally 3a values. 
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Table 2. Dust Parameters of IRS1 from the Spitzer IRS spectrum 



Model 1 Model 2 





Mass 


Temperature 


Mass 


Temperature 


Component 


(1(T 3 ^ M ) 


(K) 


(10- 3 dj M & ) 


(K) 


crystalline olivine (Mgi. 9 Feo.iSi0 4 ) 


2.3 (0.7) 


54.8 (1.6) 


2.5 (0.7) 


54.1 (1.6) 


FeO 


0.020 (0.003) 


146 (4) 


0.22 (0.08) 


89.4 (4.9) 


MgO 






0.025 (0.006) 


116 (3) 


Mgo.6Feo.4O 


0.0069 (0.0016) 


147 (6) 






amorphous silicate 


6.6 (0.3) 


57.2 (0.4) 


6.1 (0.2) 


57.6 (0.4) 


Total 


8.9 (1.0) a 




8.9 (1.0) a 





a Since the estimates of the mass of the different components are not independent, we give an upper 
limit on the error in the total mass by simply adding the errors of the individual components. 



